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Abstract 

The Goldberg- Sachs theorem is an exact result on shear-free null geodesies in a 

vacuum spacetime. It is compared and contrasted with an exact result for pressure- 

^^. free matter: shear-free flows cannot both expand and rotate. In both cases, the 

I I shear-free condition restricts the way distant matter can influence the local gravi- 

^Tj' tational field. This leads to intriguing discontinuities in the relation of the General 

Relativity solutions to Newtonian solutions in the timelike case, and of the full the- 



ory to the linearised theory in the null case. 

It is a pleasure to dedicate this paper to Josh Goldberg. 
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1 Introduction 

r> I The physical interpretation of solutions of the Einstein Field Equations in General Rel- 

c^ • ativity Theory is intimately tied in to the way timelike and null geodesies behave. The 

differential properties of families of geodesies are described by their expansion, rotation 
and shear in the timelike case [12, 15], and by the null expansion and shear in the case 
of the irrotational null geodesic congruences that underlie observations [14] . This paper 
discusses the key role of shear in physical processes as evidenced by their effect on such 
congruences, and hence the very special nature of shear-free solutions. The remarkable 
Goldberg-Sachs theorem [20] demonstrates this very special nature in the case of shear free 
null geodesies. It was proceeded by Godel's intriguing results shear-free timelike geodesies 
[19], which considered specific the case of spatially homogeneous geometries. This result 
was generalized to the inhomogeneous case of any shearfree timelike geodesies in [16]. 
Although very different in detail, the timelike and null cases are in a sense analogous 
results: they both refer to the way that shear in a congruence conveys information about 
distant matter, so shearfree congruences can only occur in restricted circumstances. That 
is what is explored in this paper. 
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The Einstein Field Equations ('EFE') take the form 

Gab = Rab — ^Rdab = I^Tab — A Qab , (1) 

showing how matter causes space time curvature by specifying the spacetime Ricci tensor 
Rab in terms of the matter stress tensor Tab- Provided the cosmo logical constant A is 
constant in time and space, the twice-contracted Bianchi identities together with (1) 
guarantee the conservation of total energy-momentum: 

{V6G'"^ = 0, VaA = 0} ^ VbT''^ = 0. (2) 

To complete the dynamical description, we must specify the matter present by providing 
suitable equations of state relating the components of Tab- 

The locally free gravitational field is given by Weyl tensor Cabcd, which is the trace- 
free part of the full curvature tensor Rabcd- It is the part of the spacetime curvature 
not directly determined pointwise by matter but rather determined by matter elsewhere 
through tidal effects and gravitational waves. The Bianchi identities 

V[e-Rafe]cd = <^ V Cabcd = —'^ [a\Rb]c — ^R9b]c\ , (3) 

are integrability conditions relating the Ricci tensor to the Weyl tensor [15, 18, 37], me- 
diating the action at a distance of the gravitational field; they take a form similar to 
Maxwell's equations for the electromagnetic field [25]. Together with suitable equations 
of state for the matter, relating the various components of Tab, equations (1), (2), (3) 
determine the dynamical evolution of the model and the matter in it. 

I will simplify things in what follows by only considering the simplest form of matter: 
a pressure-free perfect fluid, such as cold dark matter or cold baryonic matter. This moves 
with a unique 4- velocity u"" = dx'^/ds where s is proper time along the matter flow lines, 
so u^^Ua = — 1. The matter energy- momentum tensor Tab then takes the form 

Tab = pUaUb, p > 0, (4) 

where p = TabU'^u^ is the energy density. The energy-momentum conservation equations 
(2) reduce to 

p = -p Vua , iia ■-= U^VbU'' = , (5) 

so the matter, affected only by gravity and inertia, moves geodesically. This is the case 
of pure gravitation: it separates out the (non-linear) gravitational effects from the com- 
plexities of thermodynamic and fluid dynamical effects. 

2 Timelike flows 

For a given fundamental observer moving with 4-velocity tt", spacetime decomposes into 
space and time [12, 15]. The metric of the tangent spaces orthogonal to -u" is given by 

hab = 9ab + UaUb => ^hl = hi, h" a = 3, habU^ = (6) 



2 TIMELIKE FLOWS 3 

where gab is the spacetime metric. The metric hab is used to project orthogonally to m"; 
the volume element for tangent 3-spaces orthogonal to -u" is rjabc = Vatcdu'^ = Vlabc]- Here 
round brackets denote symmetrization, square brackets denote skew symmetrization, and 
angle brackets represent the Projected Symmetric Trace-Free ('PSTF') part of a tensor. 
The covariant derivative of a geodesic timelike vector field -u" may be split into irreducible 
parts as 

VfoMa = l^&hab + Vabci^" + C^ab , (^ab = Cr<ab>, (7) 

where 

Q = h'''VaUb, Cyab = KhtVi^,Ua)-'^Qhab, u'' = ^T^^^'V ^,U ^-^ (8) 

are the rate of expansion, shear, and vorticity respectively. Time and spatial derivatives 
relative to m" for a tensor T^^ are defined by: 

f\ = u'VcT\ , VcT\ = h/h^hb^Vdrf . (9) 

The Weyl tensor splits into the PSTF gravito-electric and gravito-magnetic fields 

Eab = CacbdU'^u'^ , Hab = ^VacdC^'^beU^ , (10) 

which provide a covariant description of tidal forces and gravitational radiation. 
The Einstein equation (1) with matter source term (4) and the Ricci identity 

'^la'^b]Uc = RabcdU'^ (H) 

for m" give the following evolution equations: 

e + |02 + ffafeff"^ - 2cj„a;" + i/tp - A = 0, (12) 

Cr(a6) + fQcTai + (^c{a(^bf + ^{a^b) + Eab = , (13) 

UJ{a) + lQuJa-(TabUj'' = 0. (14) 

Constraint equations are the identity VaW" = 0, the field equation 

I V'^e - Vfoff"^ - (curl w)" = 0, (15) 

where the 'curl' is (curlw)" = r^^^'^VbCUc, and an equation for the magnetic part of the 
Weyl tensor: 

ffab ^ _^{a^b) ^ ^^^^j ^yb ^^g^ 

where the 'curl' is (curlo")"^ = rf'^^'^V cO'd^ ■ Propagation equations for the Weyl tensor 
are the ii^-equation and iJ-equations: 

E<"^> -(curli/)"^ = -i/tpa"^-0E"^+ 3a<"e^^^'+ ^^"'^^"w.E/ , (17) 

H^"^^ + (curl Ef^ = - e i?'^^ + 3 a^'^c H^^^ + r]'"^^''uc H/^ , (18) 

where the 'curls' are (curlff)'^^ = t]"'^'^" V cH/^ , (curlE)"^ = r]"'^'^'' V cE/l The constraint 
equations are the (divi?) and (diviJ)-equations: 

VbE'^' = I /tVV + 3 Ub H^' + r^"^^ a,, /// , (19) 

Vbi/"^ = -Kp w" + 3 Wb E"^ - au E/ . (20) 
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These equations - the exact non-hnear equations for dust filled spacetimes - are clearly 
generalizations of Maxwell's equations for the electromagnetic field [25]. 

The case of a vacuum (empty spacetime) is a special case of these equations: just set 
p = 0. The general form of these equations for arbitrary matter fields (including pressure, 
viscosity, and heat fiux terms) is given in [18]. 

2.1 Dynamic effects 

By the vorticity conservation equation (14), the vorticity along each world line is affected 
only by the expansion and shear: for the case of pressure free matter considered here, 
vorticity cannot be created or destroyed along any world line. By the Raychaudhuri 
equation (12), the rate of expansion is increased by a cosmological constant but decreased 
by any matter present (because we have assumed the matter energy density is positive). 
This is the local Ricci effect on the fluid flow (it results from the Ricci tensor term in the 
EFE (1)), which occurs in a point-by-point manner due to the matter occurring along 
the worldline. The Weyl tensor cannot directly affect the expansion rate, but it can do 
so by inducing shear (via the shear propagation equation (13)) which then induces a 
deceleration (by (12)). This is the Weyl effect on the fluid flow; it occurs non-locally, due 
to matter at a distance from the world line. If Eat = the evolution equations along 
each world line (12) - (14) become ordinary differential equations unaffected by distant 
matter: there are no tidal effects or gravitational wave effects affecting the fluid flow, and 
each world line evolves on its own, unaffected by what is happening elsewhere. This is 
what has been called a 'silent universe' [26, 4] (for a recent review, see [42].) 

How do non-local effects occur? We can think of it as happening in two ways. First, 
tidal action at a distance is represented by the the {div E) equation (19) with source 
the spatial gradient of the energy density (e.g. scalar perturbation modes); this can 
be regarded as a vector analogue of the Newtonian Poisson equation, where by matter 
elsewhere generates an E-G.eld here. Similarly the (div if) equation (20) shows that fluid 
vorticity elsewhere generates a if-fleld here (e.g. vector perturbation modes). 

Alternatively, we can consider the evolution equations (17), (18), which together form 
a hyperbolic system. They show how gravitational radiation arises: taking the time 
derivative of the -E-equation gives a term of the form (curl H)' ; commuting the derivatives 
and substituting from the if-equation eliminates H, and results in a term in E and a 
term of the form (curl curl E), which together give the wave operator acting on E [22, 11]; 
similarly the time derivative of the if-equation gives a wave equation for H. This shows 
how matter over there can effect matter here by generating gravitational radiation which 
travels here and causes a non-zero -E-fleld here, which then affects matter here (via the 
shear equation (13). This is compatible with the constraint equation effects just discussed, 
because the constraint equations are preserved by the time evolution equations (see [40] 
corrected in [41]). As the key link between the E and H flelds in this process is via their 
curls, this suggests one can characterize the existence of gravitational radiation by the 
condition 

(curlii)"VO, (curlE)"VO (21) 

which of course requires that both E and H axe non-zero. 

However there is another way a non-zero Weyl tensor can be created where there 
was none before: this is locally via matter shear (see (17)). This emphasizes the crucial 
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importance of shear in gravitational dynamics. It is not only the link whereby information 
on surrounding inhomogeneities (given us via the electric part of the Weyl tensor) alters 
the fluid flow here, it is also a source of the electric part of the Weyl tensor. If the shear 
is zero, this link is broken, the way distant matter can influence us here is very limited: 
for by (13), only the vorticity can prevent the electric part of the Weyl tensor from being 
zero. But as we see in the next section, this does not work; distant matter is constrained 
to acting in an isotropic way around our world line. This is a very special situation. 

3 Timelike shear-free results 

When the shear is zero, the expansion is isotropic; we might expect vorticity to tend to 
generate anisotropy that would break this condition. However a non-zero Weyl tensor 
might balance this tendency. Specifically, on setting dab = in the above equations, (13) 
becomes a new constraint equation, along with the old constraint (16) determining E and 
H in terms of u: 

Eab = -ooiaCOb), H^' = -V^^u'^ . (22) 

The task now is to take take time derivatives of these constraints to see if the shear-free 
equations (obtained by setting (Tab = in all the above equations) are consistent for some 
non-trivial special cases. 

This is a major calculation: the result is not obvious. It was not initially tackled this 
way. First, in a remarkable pioneering paper presented at a world mathematics congress 
in 1950, Kurt Godel examined this question in the case of spatially homogeneous Bianchi 
IX cosmologies. He showed [19] that in this case, a shear-free universe could either expand 
or rotate, but not both; but he did not show how he had obtained that result. In 1957, 
Schiiucking derived the Godel result in detail [35]. In 1967, I used an orthonormal tetrad 
formalism to show that the restriction of spatial homogeneity was unnecessary: 

Dust Shear-Free Theorem [16]: if a dust solution of the EFE (possibly with 
a cosmological constant) is shearfree in a domain U, it cannot both expand 
and rotate in U: 

{ii" = 0, aab = 0} ^ c^e = 0. (23) 

This is an exact result, obtained by utilizing all the field equations. A covariant proof is 
given in [36]. Applying this theorem to the cosmological context, consider a shear-free 
dust-filled universe that expands. Then (23) shows Uab = 0, so from the above equations 

{aab = 0, <d>0}^Eab = 0, H''' = 0, V"e = 0. (24) 

The space-time is conformally fiat and the universe is a Friedmann-Lemaitre-Robertson- 
Walker universe [15]. These are thus the only expanding shear-free baryonic plus CDM 
cosmological solutions, provided both these components move with the same 4-velocity. 

One should note that the result (24) does not require p > 0. It is true in the vac- 
uum case (with or without a cosmological constant). Later generalizations considered 
perfect fluids rather than pressure-free matter, so acceleration of the timelike congruence 
was allowed; the result (23) remains true in all cases considered so far, indeed Collins 
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conjectured [8] that all shear-free perfect fluids obeying a barotropic equation of state 
must have either zero expansion or zero vorticity. Senovilla, Sopuerta and Szekeres [36] 
summarized results obtained towards proving this conjecture, and gave a fully covariant 
proof that shear free solutions with the acceleration vector proportional to the vorticity 
vector (including the case of vanishing shear) must be either non-expanding or non ro- 
tating. Van den Bergh [38] gave a tetrad-based approach for two particular cases require 
a special treatment, namely p + l/3p = constant, and p — l/9p = constant, as well as 
the equation of state p = (7 — l)p + constant. Van Den Bergh, Carminati, et al [39, 5] 
showed the result is generically true for shear-free perfect fluid solutions of the Einstein 
field equations where the fiuid pressure satisfies a barotropic equation of state and the 
spatial divergence of the magnetic part of the Weyl tensor is zero. 

Can one get models other than FLRW in these cases? Collins showed [7] that for 
irrotational shearfree perfect fiuids obeying a barotropic equation of state p = p{fi) and 
with nonzero acceleration, one can get spherically symmetric Wyman solutions, or models 
that are plane symmetric, and either spatially or temporally homogeneous. In all cases, 
when the space-time is sufficiently extended, the fiuid exhibits unphysical properties. 
Consequently shear-free expanding barotropic perfect fiuids must either be FLRW, or 
must be restricted to local regions where these conditions hold. Thus it turns out that 
the FRW models are the only shear-free barotropic perfect fiuid models in which the 
matter is physically reasonable globally [8]. 

Overall, these results show clearly how restrictive the shear-free result is for plausible 
fiuid models. It will of course not be true for "imperfect fiuids" with arbitrary equations 
of state: one can then just run the field equations from left to right to determine an 
unphysical form of "matter" that will give the desired result. Such calculations have no 
physical significance. One should note here that despite what one might have thought at 
first, although shear- free solutions are necessarily self-similar (they map the orthogonal 
3-spaces conformally onto each other), the converse is not necessarily true: self similar 
solutions need not be shear- free (see e.g. [6]). 

This is related to the idea of perfect fiuid as follows: a continuum description of 
matter is underlain by a kinetic theory description. Now it follows from kinetic theory for 
a collision-free fiuid that if there is non- zero shear, there will be an anisotropic stress (the 
non-zero shear will generate anisotropy in the particle distribution function which will 
then result in an anisotropic pressure [17]), hence a perfect fiuid description will not then 
apply. This will also be true if there are collisions leading to a non-zero shear viscosity, 
because of the relation 

T^ab = A aab (25) 

determining the anisotropic pressure iVab, where A is the shear viscosity coefficient [12, 15]. 
It follows that if one has a perfect fiuid, 

{7r,5 = 0,A^0}^K5 = 0}. (26) 

But kinetic theory shows that A 7^ for realistic fiuid descriptions based on kinetic theory 
with collisions [24, 17]. Hence an exact perfect fluid description of continuous matter 
implies zero shear, and the above results apply. 

The moral is that realistic matter cannot be accurately represented by a perfect fiuid 
description. The difference may be small, but is important in principle. 
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3.1 The Newtonian limit 

This shearfree result is an exact result for the full non-linear EFE: no approximations 
have been made, and is completely general: it holds for any spatially homogeneous or 
inhomogeneous model. This raises a very interesting situation as regards the Newtonian 
limit of the EFE, because this result is not true in the Newtonian case. 

The key point here is that Newtonian Gravity is not independent of General Relativity: 
it derives from General Relativity in special conditions. Specifically, it is a limiting form 
of General Relativity, valid in particular circumstances (when matter relative motion is 
at low speeds, and there are no gravi-magnetic effects or gravitational waves, which will 
be true if the magnetic part of the Weyl tensor is zero). Consequently the properties of 
Newtonian gravity should follow from those of general relativity.^ One should note here 
that to obtain Newtonian cosmological models, one has to use a description in terms of 
a gravitational potential where one allows the potential to diverge at infinity, so it is not 
strictly Newtonian theory, but rather an extension of the theory, where also the idea of 
acceleration is generalized. With these generalizations, one can obtain viable Newtonian 
cosmological models for the dynamical behavior in a matter dominated era^ (the relevant 
equations and references are given in [15]). 

The major point then is that there are shearfree solutions of the Newtonian equations 
for pressure-free matter in cosmology that both expand and rotate; specific examples have 
been given by Narlikar [27]. Consequently, the Newtonian limit is singular. Consider 
a sequence GRT{i)fj=o of relativistic shearfree dust solutions with a limiting solution 
GRT{0)cr=o that constitutes the Newtonian limit of this sequence. The latter solution 
will necessarily satisfy (23) because every solution GRT{i)„=Q in the sequence does so. 
Newtonian solutions NGT„=q that do not satisfy (23) are thus not obtainable as limits 
of any sequence of relativistic solutions GRT{j)„=q. Assuming Einstein's field equations 
represent the genuine theory of gravitational interactions in the physical Universe, this 
result tells us that not all Newtonian cosmological solutions are acceptable approximations 
to the true theory of gravity. 

An important application of this result is as follows: Narlikar has shown [27] that 
shearfree and expanding Newtonian cosmological solutions can have vorticity that spins 
up as the universe decreases in size, and hence causes a 'bounce' (the associated cen- 
trifugal forces avoid a singularity). This would be a counter-example to the cosmological 
singularity theorems of Hawking and Penrose (see [23]), if there were CRT analogues of 
these singularity-free cosmological solutions; but the shearfree result (23) shows there are 
no such GR solutions. This is a remarkable way in which the exact properties of GR dust 
solutions support the results of the Hawking-Penrose singularity theorems, obtained by 
completely different methods. 

This is a case where the Newtonian models are very misleading. The Newtonian limit 
is singular in such cases; so we need to be cautious about that hmit in other situations of 
astrophysical and cosmological interest. 



^This is in contrast to the view of some astrophysicists that Newtonian theory is the true gravitational 
theory, and GRT a set of small corrections to be made to this true theory. 
^This is not possible for a radiation dominated era. 
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3.2 The linearised case 

It is of considerable interest then whether the result (23) holds in the case of linearised 
perturbations of FLRW universe models. It has recently been shown^ that it holds in this 
case too: if a perfect fluid with equation of state p = kp in an almost FLRW universe 
is shear-free, then it must be either expansion-free or rotation-free. Thus linearization 
does not lose this property. This makes the situation even more remarkable: for then 
the linearised solutions - almost universally used to study the formation of structure by 
gravitational instability in the expanding universe, and believed to result in standard 
local Newtonian theory - also does not give the usual behaviour of Newtonian solutions in 
interesting cases. The moral of the story is that you can't believe a Newtonian solution 
unless it can be shown to be the limit of a family of GRT solutions - else it may lead you 
badly astray (as in the case of the shear-free expanding and rotating solutions found by 
Narlikar) . 

4 Null flows 

I now turn to the null case. The kinematic definitions of expansion, shear and vorticity for 
congruences of null geodesies were given by Ehlers and Sachs [14] , with a tetrad version 
being given by Newman and Penrose [28] (summaries are given in [23], [31]). This proceeds 
in parallel to the analysis for the timelike case, with two crucial differences: the geodesic 
vector field k"' = dx"" /dv is null {k'^ka = 0) rather than timelike; and the projection is into 
a two-dimensional spacelike screen space orthogonal to k"' and an observer -u", instead of 
into a 3-dimensional space as in the case of m" (so in this section, < .. > denotes trace- free 
2-dimensional orthogonal projection to k"' and u""). 

For an irrotational null geodesic congruence, the optical scalars 6 (expansion, given 
by VaA;" = 29) and aab = o'<ab> (shear) satisfy the Sachs equations 

(27) 

, .---.„ -.„ (28) 

av 

where the Ricci tensor term is $oo •= i^{p + p), determined by the matter at each point, 
and the Weyl tensor term is "^ab '■= k'^Cc<ab>dk'^, determined by matter elsewhere plus 
boundary conditions. 

These are obviously in analogy to the timelike case (12), (13) (there is no analogue to 
(14) because the vector field k"" is a gradient: ka = Va(t>, and so is irrotational). Thus one 
can again refer to Ricci focusing caused pointwise by the matter distribution $00 inside 
the beam,^ and Weyl focusing caused by the Weyl tensor term "^ab generated non-locally 
by matter outside the beam. It is the latter effect that underlies gravitational lensing and 
consequent focussing of the null geodesic rays. The full set of equations whereby these 
non-local effects take place are given in terms of the spin coefficient equations in ([30], 
pp. 248-249); these are equations (12) - (20) above expressed in terms of a null vector 



dv 


= %0, 


da,, ^ 

+ ZU aab - 


= ^ab 



•^A-M Nzoiki, R Goswami, P K S Dunsby, and G F R Ellis, in preparation. 

''it is noteworthy that the cosmological constant does not enter here: it has no direct influence on null 



focusing. 
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basis, plus equations determining the tetrad rotation coefficients and giving the tetrad 
components relative to a coordinate basis (required to get a complete set of equations). 
These equations immediately imply 

{$00 = 0,^ = 0} ^(T = (29) 

i.e. an irrotational null congruence in empty space cannot shear without either expand- 
ing or contracting. If we had included twist, the conclusion would have been altered 
to, an irrotational null congruence in empty space cannot shear without either expand- 
ing/contracting or twisting: a kind of inverse of the timelike no-shear result.^ 

4.1 The geometry of the Weyl tensor 

The Petrov classification describes the possible algebraic properties of the Weyl tensor at 
each event in a Lorentzian manifold. It was ffist given in terms of an orthonormal basis by 
Petrov [32], and is nicely described by Ehlers and Kundt [13] and by Penrose and Rindler 
([31], pp. 242-246). 

The relation to null vectors was given by Ehlers and Sachs [14], showing how the 
timelike and spacelike Weyl eigenbivectors are associated with preferred null vectors, 
called the principal null directions (PND's) of the Weyl tensor. The condition for A;" to 
be a principal null direction of the Weyl tensor is 

k[,C,],a[dkf]k'k'' = 0. (30) 

In general, there are four uniquely determined PND's at each point if the Weyl tensor is 
non-zero, but in degenerate cases two or more PND's may coincide; the degenerate PND's 
satisfy more restrictive equations than (30): k"" is a degenerate eigendirection of the Weyl 
tensor iff 

Cabc[dke]k'k' = 0. (31) 

The different possibilities lead to the six Petrov types, succinctly described using a spinor 
formalism ([29]; [31], pp. 223-226); these different algebraic types correspond to different 
physical situations. The Petrov types are, 

• Type I : four simple PNDs (generic: realistic models), 

• Type II : one double and two simple PNDs, 

• Type D : two double principal null directions (massive objects with symmetry: e.g, 
Schwarzschild), 

• Type III: one triple and one simple PND, 

• Type N : one quadruple PND: Cahcd.k'^ = (pure gravitational waves: e.g. a plane 
gravitational wave), 

• Type O: the Weyl tensor vanishes: Cabcd = (no tidal forces: e.g. a FLRW universe). 

Type I is generic, all the others are algebraically special Weyl types (they all have repeated 
PNDs) and so correspond to restricted physical situations with symmetries. 



''There is no corresponding timelike result if the cosmological constant is non-zero. 



7+^' = 

av 


= $00 


= 


= ^ab 
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4.2 The Goldberg Sachs Theorem 

As in the timelike case, for essentially the same reasons, existence of shear-free null con- 
gruences (discussed in [31], pp. 189-199)^ places strong restriction on the spacetime. We 
immediately find 

(32) 
(33) 

the first showing that no new gravitational information can enter the congruence as it 
travels from the source to the observer (only the matter encountered by the null rays 
can cause convergence) and the second shows that shear-free null geodesies are PNDs 
([31], (7.2.14) and (7.3.2)). In the vacuum case, further restrictions occur, captured in an 
important result by Joshua Goldberg and Rainer Sachs: 

Goldberg-Sachs Theorem [20]: A vacuum metric admits a shear free null 
geodesic congruence k^ if and only if k^ is a degenerate eigendirection of the 
Weyl tensor (equation (31) is true). 

This proves that a vacuum solution of the Einstein field equations will admit a shear-free 
null geodesic congruence if and only if the Weyl tensor is algebraically special. Shortly 
after the Goldberg and Sachs paper, an alternative proof was given by Newman and Pen- 
rose [28], using a tetrad or spinor formalism. A generalisation was given by Robinson 
and Schild [33], establishing a connection between algebraic degeneracy of the Weyl ten- 
sor, the existence of a null geodesic shear-free congruence, and restrictions on the Ricci 
tensor which are weaker than the requirement that there be empty space, in particular, 
they showed that the gravitational field due to any Maxwell field with shear-free rays is 
algebraically special. An even more generalized version is given by Penrose and Rindler 
([31], pp. 195-198). 

The theorem is useful in searching for algebraically special vacuum solutions, which is 
very helpful because almost all solutions we can write down in exact form are algebraically 
special, corresponding to restricted matter distributions and boundary conditions; exam- 
ples are the Kerr solution and plane gravitational waves. 

4.3 The News function 

Given the discussion above, one might expect that shear relates to the way information is 
conveyed along bundles of null geodesies, with restrictions on that information when the 
shear is zero. This expectation seems to be borne out by the analysis of axisymmetric 
vacuum spacetimes by Bondi, van den Berg and Metzner [3] , conveniently summarised in 
the book by D'Inverno [10]. On using a null coordinate system where the past null cones 
of the central observer are given by {u = const} where u is the retarded time and r is a 
measure of distance down the past light cone, the metric is the Bondi metric 

ds^ = -{-e^f -r^e^''U^)du^ + 2e^^dudr + 2Ur^e^^dude-r\e^"'de^ + e-^"' sin^ ed(f)^) (34) 
r 



^Shcar free null congruences play an important role in twister geometry ([31], Chapter 7). 
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where V = V{u,r,9), U = U(u,r,6), [5 = V{u,r,9), 7 = 7(-u,r, ^). Then solving the 
vacuum EFE asymptotically, V = r — 2M + 0{r^^), /3 = —n^/Ar"^ + 0{r^^)^ Aq^^ = 
2Mn — do + d cot 6, and 

7=^^^^ + ^ + 0(r-). (35) 

The mass of the system as measured at infinity is the Bondi Mass 

m{u) = - f M{u, 6) sin OdO (36) 

The shear of the radially outgoing null rays is 

a = n{u,e)/r'^. (37) 

The initial data is 7(m, r, 9) on an initial value null surface, plus n^u{u, 9) which determines 
the evolution of the source and so is called the News Function^ determined by the first 
time derivative of the shear. Finally 

mo = -- r(nofsm9d9 (38) 

2 Jo 

which is non-positive and so shows that there is mass loss if and only if there is news. 
The Weyl tensor up to order 1/r has the non-zero outgoing radiation component 

^0 = -n^uu/r = -r a^uu (39) 

so Bondi's formula says that the Bondi mass of the system decreases if and only if there 
is outgoing radiation, and so if and only if the second time derivative of the shear is 
non-zero. As stated in the abstract, "/t is shown that the flow of information to infinity 
is controlled by a single function of two variables called the news function. Together with 
initial conditions specified on a light cone, this function fully defines the behaviour of the 

system The principal result of the paper is that the mass of a system is constant if and 

only if there is no news; if there is news, the mass decreases monotonically so long as it 
continues." This confirms the key role played by shear in conveying news, and hence in 
the mass loss carried out by outgoing gravitational radiation. The result remains true if 
the assumption of axisymmetry is dropped [34] . 

Thus we have confirmation of the key concept of this paper: the shear of geodesies, 
generated by the Weyl tensor, conveys information about the gravitational field due to 
gravitating bodies. If the shear is zero, this information is very limited and the space time 
dynamics is highly constrained. 

4.4 The discontinuous limit 

However there is an intriguing context where this picture seems to fail. The outgoing grav- 
itational waves discussed in the previous section should at large distances asymptotically 
become plane gravitational waves. 

The pp wave geometries [1] are given in terms of null coordinates by 

ds^ = —2dudv + hij{x, u)dx'^dx^ (40) 
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with a covariantly constant null vector field C,"' = (d/du)"" which is consequently a Killing 
vector whose expansion and shear vanish. Hence they are subset of null shear-free solu- 
tions; the Weyl tensor is type N and this vector field is a 4-fold degenerate PND (thus 
satisfying the Goldberg-Sachs theorem). Plane gravitational waves [2] are a special class 
of vacuum pp-wave where 

ds"^ = -2dudv + [a{u){x'^ - y"^) + 2b{u)xy]du'^ + dx^ + dy"^ (41) 

Here, a{u) and h{u) can be any smooth functions; they specify the amplitude of the two 
polarization modes of gravitational radiation. The waves can convey arbitrary messages 
by the time variation of these modes along successive light cones. 

The issue then is the following: 

1: The Bondi news function for asymptotically flat vacuum metrics is the time deriva- 
tive of the shear. Hence, no shear implies no news. 

2: The BMS metrics (34) should become plane gravitational waves asymptotically at 
infinity. 

3: The pnd vector field of a plane gravitational wave is necessarily an exactly shear- free 
geodesic congruence (Goldberg-Sachs). 

4: But plane fronted gravitational waves can freely carry news: they have two free 
polarisation functions, even though their shear is exactly zero. 

5: Hence the Bondi metrics should tend to a state where there is no shear at infinity, 
hence no news, yet effective news transfer is possible in the exact limiting spacetime. 

How does this all fit together? Presumably it is a question of relative orders of decay 
near infinity: but it is not obvious how it is coherent! One might have thought that zero 
shear meant zero news: but this is not the case! 

There is of course a shear associated with plane waves, which is essential for the 
existence of the waves. Null geodesies which intersect the null hyperplane histories of 
the plane waves must have shear (if they don't then the plane waves don't exist). Other 
manifestations of this phenomenon are Penrose's observation that an impulsive plane 
gravitational wave acts as an astigmatic lens and also that colliding plane waves generate 
shear after collision. But the paradox remains: the PND congruence is shear free, and 
indeed as a consequence the geometry is very limited (it has a large symmetry group) but 
still can convey arbitrary information. 

4.5 Linearised gravity 

There is no Newtonian analogue of the Goldberg-Sachs theorem, as there are no equivalent 
concepts there. But one can consider the linearized version of the result. It has been shown 
by Dain and Moreschi [9] that a corresponding theorem will not hold in linearized gravity, 
that is, given a solution of the linearised Einstein field equations admitting a shear-free 
null congruence, then this solution need not be algebraically special. 

This is a warning of the perils of using linearised results for a non-linear theory: some 
key results may not be valid in the linearised case, even though they are an exact result 
of the non-linear theory. Hence as in the case of the Newtonian limit of the timelike case, 
here we can consider a sequence of exact shearfree solutions that tend to a linearised shear 
free solution: the limiting properties of the exact solutions differ from the properties of 
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the linearised solution. Hence for example special solutions of the exact and linearised 
equations may have differen properties. 

5 Conclusion 

Shear of geodesic curves plays a crucial role in conveying information regarding the state 
of matter in a region, underlying structure formation in the timelike case and gravitational 
lensing in the null case. 

I have compared the timelike and null cases of shear-free geodesies. In each case 
the allowed exact solutions are strongly restricted, albeit in rather different ways. The 
argument makes clear that the two contexts are related in the following way: one would 
not expect gravitational radiation to be emitted by a shear-free dust flow, and that is 
indeed the case (combine (24) and (21)). We have also seen that there are intriguing 
limiting questions in each case: shear-free solutions don't have the same implications in 
Newtonian theory (timelike case) and linearised gravity (the null case). Pursuing these 
issues might be interesting. 

One further line of of investigation that may also be interesting is to pursue the sta- 
bility of the Goldberg-Sachs result in the following sense: an "almost Birkhoff theorem" 
has recently been proven [21]. This shows that if the conditions for Birkhoff 's theorem 
are almost true, then Birkhoff's theorem will be approximately satisfied. An "Almost 
Goldberg-Sachs theorem" would similarly show that if the conditions for the Goldberg- 
Sachs theorem are almost true, then Goldberg-Sachs theorem will be approximately sat- 
isfied. If this were not true, that might throw light on the failure of the asymptotic limit 
discussed in the previous section. 

It is a pleasure to dedicate this paper to Josh Goldberg, who has done so much for GRG 
studies. I thank Peter Hogan for helpful comments, and Roy Maartens for suggestions 
that have greatly improved this paper. 
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